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We introduce an effective scalar field theory to describe the ''He phase diagram, which can be 
^ ■ considered as a generalization of the XY model which gives the usual A-transition. This theory 

O ' results from a Ginzburg-Landau Hamiltonian with higher order derivatives, which allow to produce 

transitions between the superfluid, normal liquid and solid phases of *He. Mean field and Monte 
Carlo analyses suggest that this model is able to reproduce the main qualitative features of *He 
phase transitions. 
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I. INTRODUCTION 



, Effective field theory models are widely used to describe phase transitions in condensed matter systems. The 
^ ■ simplest example is that of the Ginzburg-Landau theory T], in which a "Hamihonian" Hohl't'i] depending on certain 
O field variables (pi and defined on the sites i of a lattice or at points x, 4){x), in its continuum formulation, may account 
for the description of first or second-order phase transitions. 

The effective theory is not intended to be a microscopic theory at aU. Instead, it makes use of the fact that critical 
phenomena are divided into universality classes which are determined by a few basic properties of the system only, 
such as the dimensionality of the space, the range of interactions, the number of components and the symmetry of 
the order parameter. The renormalization-group theory predicts that, within a given universality class, the critical 
exponents and the scaling functions are the same for all systems so that we can make use of the corresponding simpler 
effective theory to calculate such quantities. 
■ In this context, the superfluid transition of ^He, occurring along the A-line T\{P), where P is the pressure and 
T the temperature, belongs to the three-dimensional XY universality class 0- Its order parameter is related to the 
; complex quantum amplitude of helium atoms so that the 0(2) field theory may serve as an effective description 
' of this transition. In fact very good agreements between the critical exponents and scaling functions of this model 
^ , and the experimental measurements of the A-transition of ''He are found 0| . 

^ Fig.Qshows the ''He phase diagram in the (T, P) plane. The 0(2) field theory is intended to describe this system 

in the vicinity of the A-transition, and strictly speaking only near at zero pressure, because the model presents a 
temperature-driven transition only. The whole A-line is however expected to belong to the same universality class for 
Q ] P 7^ 0. The transition lines between the liquid phases and the solid phase are experimentally observed to be of first 
Q . order, with a finite entropy difference between the phases and the presence of a latent heat j^]- 

In Fig.^ the solid, superfluid and normal liquid phases meet at a single point, suggesting the presence of a Lifshitz 
point ^] in the ^He phase diagram. In order to examine the critical behavior around this point, it would be very 
' convenient to have an effective theory of this system containing all the three phases which meet there. This model 
^ ' should then be able to account for the transitions between the solid and the superfluid phases, solid and normal liquid 
- - < phases, and superfluid and normal liquid phases. 

We address the problem of finding such an effective theory in the present work. We take as an starting point the 
0(2) field theory, which describes the transition between the superfluid and normal liquid phases. In this theory, the 
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FIG. 1: Experimental pressure-temperature phase diagram for *He in the low temperature region. Two liquid phases, He I and 
superfluid He II, and a solid phase can be distinguished. 



superfluid phase corresponds to a ferromagnetic phase, in which there exists long-range order, while such an order is 
absent in the normal liquid phase, which corresponds to an unordered or paramagnetic phase. In order to generalize 
the 0(2) field theory to a model including also a solid phase we should extend the regime of applicability of the 
effective theory towards shorter length scales below the lattice spacing of the solid. 

What is the salient feature of the interactions at distances comparable to this lattice spacing? The periodic ordering 
in the ground state underlines the importance of the distance dependence of the forces acting in this scale regime. What 
kind of terms of the effective theory can reproduce an important distance dependence? The ultralocal polynomials of 
the field, without space-time derivatives, represent momentum (i.e. distance) independent interaction strength. When 
such a vertex is inserted into a graph then the distance dependence of the propagators can induce distance dependent 
vertex functions and therefore forces. But this is an indirect mechanism, the distance dependence actually originates 
from the propagators, the usual O (p^) part of the action. When the vertices carry explicit momentum dependence 
then the resulting forces should display more pronounced distance dependence. Therefore one suspects that terms 
with higher order derivatives might be the key to reproduce the solid phase. 

The simplest vertex with higher order derivatives is quadratic in the field variable, meaning the modification of the 
free dispersion relation. If the dispersion relation turns out to be negative at a certain momentum, then the elementary 
excitation with such a momentum starts to condense in the vacuum. This condensation will be stopped and the system 
will be stabilized by the repulsive O {(f>'^) interaction. Starting from the normal phase, the appropriately chosen mass 
term or the coefficient of the higher order derivative term will take us into the superfluid or the solid phase which will 
be called modulated phase. Note that the effective model may also contain an antiferromagnetic phase, in which the 
period length is the shortest possible scale of the theory. Such a phase will be however non-physical, since effective 
theories are supposed to give sensible physical descriptions only at energies lower than their energy cutoffs. 

The universality argument of the normal-superfluid transition can be extended towards the solid phase at least in 
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the leading order of the perturbation expansion. In fact, the perturbative, one-loop renormalizability at the Lifshitz- 
point |7| provides us the universality at this order. It remains to be seen if the universality can be established 
nonperturbatively. 

The heuristic reasoning above applies to any system with modulated ground state, such as the usual solid state 
crystal, Wigner-lattice of dilute electron gas and charge density wave state. But there is another consideration which 
makes the argument plausible in the context of "^He, namely the existence of rotons. The higher order derivative 
terms of the effective action are present in either phases, their strength varies only when the phase boundaries are 
crossed. As we will see, the existence of rotons, the enhancement of elementary excitations at a given momentum 
which corresponds to a local minimum in the dispersion relation, is naturally reproduced by means of an action 
whose quadratic part in the field contains higher order derivatives. In fact, the rotons should correspond to the local 
minimum of the dispersion relation which becomes the absolute one as the superfluid-solid transition is crossed. 

It is interesting to notice the formal similarity between the dynamics which drive the normal-to-superfluid and 
the normal-to-solid transitions. The former is spontaneous symmetry breaking where the potential energy reaches 
its minima at nontrivial, i.e. nonsymmetrical values of the field, at a nontrivial scale in the internal space. The 
symmetry spontaneously broken is an internal one and the dynamics is modified mainly in the IR domain. Similar 
phenomenon may take place in the external space, which leads to the modulated phase. If the dispersion relation 
reaches its minimum at nonvanishing momentum then particles with such momentum condense and the vacuum 
becomes modulated. This mechanism is driven by the derivative terms and modifies the dynamics mainly at length 
scales comparable with the inverse momentum of the particles condensing. The symmetry broken dynamically is an 
external one, rotations and translations. 

The organization of our paper is the following. In Section II the effective theory for helium-4 is defined both in 
mathematical and physical terms. Section III is devoted to the mean field solution of the model. In Section IV 
we compare the mean field predictions with numerical Monte Carlo simulations of the system. We analyze its 
phase diagram and investigate the order of the transitions between phases. Finally, the conclusions are presented in 
Section V. 



II. EFFECTIVE THEORY 



Our model is an extension of the Ginzburg-Landau model in three dimensions for a complex scalar order parameter 
by adding higher order derivatives to the action. 
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A. 



(1) 



where j = 1, 2, there is an implicit summation over repeated indices, and the kinetic energy contains the function 

A:(z) = 1 + C2Z + C42^ (2) 

The theory is regularized on a lattice of spacing a (the ultraviolet cutoff, which will be taken to be smaller than the 
period of the solid phase) , so that the integral becomes J^x ■ ^^'^ ^^^^ ^^"^ dimensionless lattice field variable 

4'x = \/0'4>{x), and gets 



E 



(3) 



with A4>x = J2 i4'x+(i + (t>x-fi -24>x), ^ mla? , A = aAc 
Expanding A/C(A) in S'[(/)], we obtain 



S[4>] = ~KlSl - K2S2 - K3S3 - K4S4 - K5S5 - KeSe + K'^(t)j,x4'j,x + ^ E i^3,x4>3,x) 



(4) 
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with the following definitions for the different terms and coefficients: 

Si ^'^4'j,x(t>j,x+ii, Ki = 1 - 12c2 + 123c4 , (5) 

5'2 = ^ 0j,2;0i,a;+2A 7 = C2 - 18 C4 , (6) 

53 = ^ (pj^xf^j^x+Zjl, K3 = C4 , (7) 

5*4 = ^ (0j\a:+A+i> + 0J,:r-A+i>) ' ^^4 = 2 C2 - 36 C4 , (8) 
x.^<iv 

S5 — {(t>j,x^2fl+i> + 4^j,x-2fL+i> + </'j,2;+2i>-f /i + 0j,a;-2i>+/i) , ^^5 = 3 C4 , (9) 

S'e = X! '^^J^^^ ('?^j,x+i+2+3 + '^j,a;+i-2+3 + 0j,:E-i+2+3 + ^j.a^-l-S+s) ' Kg = 6 C4 , (10) 



K= — + 3-21c2 + 162c4. (11) 



Physical interpretation 



This effective model presents a number of parameters: m^, A, C2, and C4. and A offer the possibility to have a 
broken symmetry, as in the standard Ginzburg-Landau Hamiltonian. Indeed, the polynomial 

Hiv) = ^rQifi^ + ^"0^*, uo>0 (12) 

has two minima for tq < 0. The coefficient uq must be positive so that lim;^^±oo H{ip) = +00, which guarantees the 
stability of the minimum of H{ip). 

Therefore, < and > will give the broken and the paramagnetic phases, respectively, in the tree-level 
approximation. The parameter has to vanish linearly at the critical temperature, cx (T — Tc). 

The other two parameters, C2 and C4, appearing in the higher order derivative part of the action, will be responsible 
for the emergence of the solid phase. Since making C2,C4 in Eq. |^ gives the usual 0(2) field theory, which 
describes the ^He transition at P = 0, we will associate the pressure to a linear combination of these two parameters, 
to be determined afterwards. 

III. MEAN FIELD SOLUTION 

The first, simplest step in determining the phase structure of the model is the tree-level solution which produces 
an inhomogeneous mean field in our case. Since this model is translational and rotational invariant in space we shall 
look for a mean field vacuum of the form 



n,x 



/cos (i^'^x'^) 



(13) 



where the amplitude (f> and the numbers K^^, fj, = 1,2,3, serve as variational parameters to minimize the action. 

The action density, s = S/L'^, for an homogeneous, ferromagnetic vacuum {Rf^ = V/x) on a lattice L"^ is obtained 
by minimizing 

SFM = + (14) 



-3m4/2A 







<0, 
> 0. 



(15) 
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In order to study the general case K^^ ^ 0, we shall need the eigenvector of the lattice box operator, 

A(j),,, = -P'^c^j,,, (16) 

where 

P2=45:sin^(^). (17) 

One finds 



-A/C(A)0j-, =A^2^,-^^ = P''{l-P''c2+P^Ci). (18) 

Then the mean field action is 

s = \{m^+M^)cl? + ^<t>\ (19) 



with 



minimum 



__3K + Af)! ,2 __6(m2 + Af) 

where + = m? + P'^{1 — P^C2 + -P^C4). This solution is vahd only if + M.^ < 0. The extrema of to^ + 
are reached at 



^2 _ C2 ± \/cl- 3C4 



^ ^ 1 (21) 



for C4 7^ 0, and 



1 

2^2 



Po^T^ (22) 



for C4 = 0. However, P^ is a local minimum of Eq. (|20|l . while is a local maximum. 

The local minimum will be global or not depending on the specific values of C2 and C4. Moreover, could 
be larger than 12, which is the maximum allowed value for P^ from its definition Eq. (|17|) . or it could even be an 
imaginary number. It is therefore necessary to carry out a careful analysis of the absolute minima of Eq. H2U|I . which 
will give us the different mean field vacua in the (c2, C4) plane. 

A. Ordered mean field phases in tlie (02,04) plane 

The conditions that will determine the global minimum of Eq. ()20|l for given (c2, C4) are: 

Existence of Pi C4 < — (23) 



P7 < 12 C4 > max 



3 

"C2 1 



(-1 + 24C2) 
36 432^ ' 



(24) 



M^{P^ = 12) < ^ < Y^(-l + 12c2) (25) 
X2(p2 = p2) < ^ C4 < ^. (26) 

In terms of these conditions, and supposing a value of such that + A^'^ < 0, the phases (values of K^) which 
minimize Eq. (|20|) are (see Fig. |21for an specific example of every case): 

1. If condition H23|) does not hold, then Eq. H2U|) is minimized by = V/z: ferromagnetic vacuum. 

2. If condition 1)23(1 holds, but condition (|24|1 does not hold, then 
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FIG. 2: as a function of for different values of (02,04) corresponding to the cases signalled in the text (from left to 
right, top to bottom): 1, 2 (a), 2 (b), 3 (a) and 3 (b). 



(a) If condition (|25|l does not hold, then the minimum is again ~ V/i; ferromagnetic vacuum. 

(b) If condition (|25|) holds, then the global minimum of Eq. (|20|) is at ~ 12 or K^^ = tt V/i: antiferro- 
magnetic vacuum. 

3. If both condition ^2[\\ and condition H24(l hold, then 

(a) If condition (|26|l does not hold, then the minimum is the ferromagnetic vacuum. 

(b) If condition (|26|) holds, then the minimum of Eq. (|2(J|I is at the vector such that — P^: we will call 
this a modulated phase. 

Note that the plot in Fig. El case 3(a), corresponds to the dispersion relation of phonons and rotons (Landau 
spectrum Q): we are in the ferromagnetic phase (superfluid phase) and there are two kind of excitations, at zero 
and different from zero momenta. At low energies, the curve is a straight line, corresponding to a phonon dispersion 
relation, while at higher energies, the spectrum deviates from a straight line, passing first through a maximum and 
then a minimum. The excitations with energies near this minimum are called rotons. The existence of the finite 
energy gap for rotons is crucial for the superfluidity in He II. The system enters into the solid phase when this gap 
tends to zero. The shape of the Landau spectrum has been confirmed by neutron-scattering experiments carried out 
in several different laboratories |^] . 

From Eqs. H23|) - H2t)() and the previous discussion, we can obtain the range of values of (c2, C4) for which the mean field 
vacuum is ferromagnetic (FM), antiferromagnetic (AF) or modulated (MOD). The result is summarized in Tabled 
and the phase diagram which results in the (c2, C4) plane is plotted in Fig. 13 

Let us now consider the order of the FM-MOD and MOD-AF transitions in this mean field approach. To do so, 
we will study the variation of the minimum of the mean field action Eq. H2U|) along a line C2 = const. Let us take 
the line C2 = 0.5 which is plotted vertically in Fig. O The value of Smin along this line depends on A4^, which in 
turn depends on the value of P^ = P^ which minimizes it. This value is P^^ — 12, i-+ and 0, respectively, for the 
antiferromagnetic, modulated and ferromagnetic phases, and is plotted in Fig. ^ along the C2 = 0.5 line. From this 
figure one can see that P^ has a discontinuity when passing from the MOD to the FM phase. This seems to suggest 
that Smin is also discontinuous and that the transition FM-MOD will be first-order. However, we will now see that 
this is not the case. 

The value of Smin at the FM phase is given by (see Eq. ((JOl, with ~ 0): 

FM _ (2'7) 
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FM 


V C4 


— < C2 < 

12 - 


1 
6 


AF 
FM 


ifc4 < ^(-1 + 12C2) 
ifC4 > ^^^(-1 + 1202) 



- < C2 < OO 

6 



AF ifc4 < — (-1 + 24C2) 

1 

MOD if — -(-1 + 24C2) < C4 < 

432 4 

FM if C4 > — 

4 



TABLE I: Mean field phases for given (c2,C4) values, supposed that + A^^ < 0. 
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FIG. 3: Mean field phase diagram in the (c2,C4) plane, supposed that rr? + M'' < at every point of the plane. The black 
zone represents the AF phase, the grey zone is the MOD phase, and the white region, the FM phase. 

This value is independent of C2 and C4. The value of Smin at the MOD phase is given by Eq. H2U|) . with M?' = A4'^{P^). 
This gives 



MOD _ (2 4 + 2(4- 3c4f/^ - 9 C2C4 - 27 elm 
" 486 c|A 



(28) 



However, when we evaluate this expression along the curve separating the FM and MOD phases, C4 = C2/4, we obtain 



'(C2,C4 = C2/4) = - 



3m 
~2X 



FM 



(29) 



and there is no discontinuity in Smin between the two phases. Therefore, we conclude that the transition FM-MOD 
is a second-order transition. 

In fact we could already have arrived to this conclusion just by a careful inspection of Fig.|21 The transition between 
the FM and MOD phases in Fig. O corresponds to the transition between cases 3 (a) and 3 (b) in Fig. [3 In case 3 (a), 
M'^{P1) > and the minimum of the action is at = (A^^ = 0). In case 3 (b), A4^{P^) < and the minimum of 
the action is at P^. At the transition, = and there is no discontinuity in the action, but there is in its derivative 
with respect to the parameter of variation. The transition is second-order. 

With respect to the AF-MOD transition, taking place at the curve C4 = (—1 + 24c2)/432, it is evident from Fig.^ 
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FIG. 4: Values of which minimize at C2 = 0.5 as a function of C4. They define the different vacua of the mean field 
action, as explained in the text. 



that Smin should be continuous since P^^ is continuous at this transition. In fact, again from Eq. (|20|l . 



Ap 3(12- 144 C2 + 1728 C4 



m 



2\2 



2A 



(30) 



and at the transition line 



AP ( {-l + 2Ac2) \ mod( {-l + 2Ac2) \ 3(8- 48 C2 + m^)2 

C4 = — j = [C2,C, = — j = . (31) 

This transition corresponds to that between the cases 2 (b) and 3 (b) in Fig. [3 In the first case, > 12 and then 

the minimum is at Pf^ — 12, A^^(P^j = 12) < 0. When we approach the modulated phase then 12. At the 

transition = 12 and therefore changes in a continuous way between the two phases. In this case the derivative 
of is also continuous. 

Finally, the transition between the AF and the FM phases in the second region of Tabled {^j^ < C2 < ^), turns out 
to be also continuous, in spite of the fact that Pj^^ jumps from 12 to at the line C4 = (— 1 + 12c2)/144 since 

AF ( _(-l + 12c2)\ 3m4 _ pj^^ 



sZn (^C2, C4 = ^ j = = s^t^ (32) 

This can again be understood as the transition between cases 2 (a) and 2 (b) in Fig. 12 when we pass from Ai^ — 
at the minimum of the action in the first case, to < in the second case. The change in M."^ is then continuous, 
but its derivative is discontinuous. 

In conclusion, all the transition lines in Fig. O turn out to be continuous transitions. At the FM-AF and FM-MOD 
transitions, when there is a jump in P^, the derivative of the action with respect to the parameter of variation in the 
(c2,C4) plane is discontinuous, while this derivative is continuous at the AF-MOD transition. 



B. Complete mean field phase diagram 

In the previous analysis we have made the assumption that was such that m? + M.^ < 0, which guarantees 
that we are in the broken phase, i.e. (p^in > (see Eq. (I20|l '). If this condition does not hold then the minimum of 
the action (|19|l is at = and we are in the paramagnetic (PM) phase. This happens at a different value of to^ 
depending on which ordered phase is considered: 

{FM phase < 

AF phase < -M^{12) (33) 
MOD phase < -M^{Pl) 

In the three dimensional plane (m^, C2, C4) we will have then four different phases: PM, AF, FM, and MOD. As we 
argued in sections HI and ITU we are interested in the physical region containing only FM, PM and MOD phases, which 
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FIG. 5: Mean field phase diagram in the (m^,C2) plane, where at every point C4 = 0.25c2. Here the phases are: PM (white 
region), FM (grey region), and MOD (black region). 

will represent the superfluid, normal liquid and solid phases, respectively, as a function of two parameters: m^, which 
would correspond to the temperature T, and a combination of C2 and C4, which would correspond to the pressure 
P. Taking C4 proportional to C2 with an adequate slope gives us a section in the (02,04^) plane represented in Fig.|3| 
where the qualitatively correct phase diagram appears: a FM phase followed by a MOD phase. The exact value of 
the slope is arbitrary, and we make the choice C4 = 0.25c2 so that the transition line between the MOD and FM 
phases is located at C2 = 1. The complete phase diagram we get in this way is shown in Fig. [S] in the (m^, C2) plane. 
Note that this physical region may not be stable under renormalization group (RG) transformations. Under these 
transformations, which keep the physics fixed (unmodified long distance behaviour), the period length will appear 
shorter when expressed in lattice spacing units. This means that they will in general connect the MOD phase with an 
antiferromagnetic phase, but with a quite complicated blocked action containing much more parameters than simply 
C2 and C4. Therefore the non-physical AF phase which our simple effective model contains will not be connected by 
RG transformations with the physical region and we can safely discard this phase and its vicinity from our analysis. 

It is easy to see that the transitions from an ordered phase (FM or MOD) to the paramagnetic phase are also 
continuous. The action inside the PM phase is in the mean field approximation while in the FM and the MOD 
phases it is given by Eq. 119|) . However the transition lines to the PM phase are just the regions of the parameter 
space where + ^ and consequently the action is vanishing. 



The preceding mean field study is qualitative only due to the absence of fluctuations. In order to estimate this 
error and to have a more reliable phase structure a numerical simulation of the model was performed as well. We 
considered a three dimensional cubic lattice of side L with periodic boundary conditions and studied the model 
Eq. Q as a function of the two parameters and C2, taking A = 0.1 and C4 = 0.25c2 as fixed values, in lattices 
L = 6,8,12,16,24 with a Monte Carlo simulation. The algorithm was a standard Metropolis, and the errors were 
computed with a jack-knife method. In the largest simulations, we performed up to 50 million of full-lattice sweeps 
(measuring every 5 sweeps). We checked that the autocorrelation times were small with respect to the number of 
measurements performed for every lattice size. 



IV. MONTE CARLO SIMULATION OF THE MODEL 
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A. Observables and the phase diagram 



In order to identify the possible inhomogeneous condensates we write the field variable 4>x in terms of its Fourier 
transform, 

0. =^e'^"0/^, ^^, = l^e-^-</,,, (34) 



K 



where Kx = K^x'^\ and there is summation over the repeated index fi. The possible values of K on each lattice 
direction arc (27r/i)n, with n — 0, . . . , L — 1. The magnetization corresponding to the wave-vector K is defined as 
the magnitude of the Fourier coefficient, 



Mk^\/UkP. (35) 

Our mean field solutions are pure Fourier modes, i.e. such that H'/'xlP — ^ ii^o for a certain which 
minimizes the action, and we introduce the corresponding momentum square as 

P2=4^sin^(^). (36) 

is the only combination of the components of the wave-vector K^^ appearing in the mean field solution. But this 
means that all different K configurations with the same value of (K) are degenerated in energy. The degree of this 
degeneration naturally may depend on the value of P^. The surface in if-space which corresponds to a given value 
of P^ is depicted in Fig. Elfor three different values P^ — 1, 3, 6. 

The configurations are not single Fourier modes of the mean field type in the Monte Carlo simulation. A possible 
generalization of P^ is the average momentum square over all modes, 

^^4^^||^^|psin^f^). (37) 



K 



Another useful observable is the energy or in general any of the eight terms (or any linear combination of them) 
appearing in Eq. Considering the energy as a function defined on the plane (m^, C2), the direction of the fastest 
change at a transition will be the one orthogonal to the transition line. Since the FM-MOD and FM-PM lines are 
almost horizontal and vertical in the parameter space, respectively, we considered the coefficients of C2 and in 
S. Both linear combinations proved also to give a good signal for the MOD-PM transition. They can be read from 
Eqs. (El) and fK|-fTT)l: 

Sni = ■^^<t>3,x<j)j,x, (38) 

X 

Sc = 56.25 Si - 10.5 S2 + 0.75 S3 - 42 S4 + 9 S5 + 6 Se - 39 . (39) 
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FIG. 7: Monte Carlo phase diagram for an L = 16 lattice (continuous line and dots) versus mean field phase diagram (dashed 
line). Error bars which are not shown are smaller than the size of the points. 

These energies are also the appropriate ones to extrapolate the mean values of observables in an interval around a 
certain simulation point by means of a Ferrenberg-Swendsen reweighting method [13. The observables Eqs. ^ 
and H35|l . together with the appropriate energy terms were measured in Monte Carlo simulations. 

We performed two different kinds of simulations. First we swept the whole parameter space by fixing the value of 
one of the two parameters (c2 or m^), varying dynamically the other one in small steps (typically 10"'^) and measuring 
observables after a number of iterations (around 5000 Monte Carlo steps). This procedure allows us to locate the 
transition point by means of the rapid changes experienced in the different observables. We call this kind of simulation 
an hysteresis owing to the typical signs of metastability observed when crossing a first-order transition (see Figs. El 
and 1101 below). Once the phase diagram had been outlined we performed better statistics simulations at fixed values 
of the parameters at the transition lines to get a deeper insight into the properties of these transitions. 

Figure 13 shows the observed phase transitions in an L = 16 simulation, together with the mean field lines. One can 
see that the Monte Carlo results agree well with the mean field phase diagram. It is interesting to see that near the 
point where all phases meet, the slopes of the different transitions are similar to those of the helium phase diagram, 
see Fig. ^ after making the correspondence of the MOD, FM and PM phases to the solid, liquid He II and liquid 
He I phases, respectively (but note that the model does not pretend to give the exact slopes of the transition lines). 
This effective model seems therefore to give a good description of the helium phase diagram around the Lifshitz point 
when C2 and are considered to be proportional to the pressure and the relative temperature, respectively. In the 
following section we study the nature of the phase transitions by using a Monte Carlo simulation and confirm that, 
in contrast to what we saw in the mean field analysis, they are of the same type as those of helium. 



The FM-PM transition is clearly second order, as we will show by calculating its critical exponents. This behavior 
is expected in order to reproduce the A-line that separates the liquid He I and He II states in the phase diagram of 



Here the order parameter is the standard magnetization Mq and one can perform the usual finite size scaling 
analysis in order to extract critical exponents and critical temperatures. A good choice for the scaling variable is the 
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0.158 



-0.156 



FIG. 8: C/i for C2 = 0.5 and lattice sizes L = 8, 12, 16, 24. 



correlation length. In a finite size lattice it can be defined using a second moment metliod [ll| as 



1/2 



4sin^(7r/i) 



(40) 



where Mq and are defined in Eq. (|35|l . and if„j = {2tt/L, 0, 0) is the minimum wave vector compatible with the 

periodic boundary conditions. 

For an operator O that diverges as where t is the reduced temperature the mean value at a temperature T 

in a lattice of size L can be written in the critical region by means of the finite-size scaling ansatz 



OiL, T) = L^olu (Fo(e(i, T)IV) + 0{L-^)) , 



(41) 



where Fq is a smooth scaling function and uj is the universal leading correction-to-scaling exponent. In order to 
eliminate the unknown Fq function we use the method of quotients ^3 where one studies the behavior of the 
operator of interest in two lattice sizes, L and rL, 



Qo^O{rL,t)/0{L,t) 



(42) 



and one chooses a value of the reduced temperature t such that the correlation-length in units of the lattice size is 
the same in both lattices. This temperature can be considered as the apparent transition point for the size L. One 
obtains easily 



Qo\ 



0{L- 



(43) 



We used the quotient method for pairs of lattices of sizes L and 2L and determined the values of the parameters 
(m^,C2) where the f/L curves cut each other. This is shown in Fig. |Slfor C2 — 0.5 and i = 8, 12, 16, 24. The critical 



exponent v (^ • 



') was measured by using 9rn2^ {dm^S, 



as the observable O of the quotient 



method and the apparent exponents together with the transition points are shown in Table Ull The calculation of the 
scaling corrections and the exact extrapolation of the critical point and exponents have not been carried out since our 
intention was just to check that the transition belongs to the university class of the XY model in 3 dimensions (which 
corresponds to C2 =0). More precise calculations of critical exponents would require the use of update algorithms 
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L 


mtL 




6 


-0.1562(6) 


0.663(10) 


8 


-0.1613(4) 


0.656(12) 


12 


-0.1629(5) 


0.678(17) 



TABLE II: determined by the intersection of the correlation lengths measured in two lattices of sizes L and 2L and the 
apparent critical exponent v obtained from the quotient method applied to the same (L, 2L) pairs. 

with smaller autocorrelation times in the vicinity of a continuous transition than the standard Metropolis, such as a 
single-cluster algorithm (for the problematics of the application of these algorithms to the standard Ginzburg-Landau 
model see e.g. Rcf. 15]). For us it is enough to confirm that the values of the exponent v reported in Table UTI turn 
out to be fully compatible (as expected) with that of the XY model in Sd 0.67155(27). 

2. FM-MOD transition 

Hysteresis-type simulations along fixed rn? < lines show metastability signs, indicating a first-order character of 
the FM-MOD transition at C2 ^ 1. This is shown in the top part of Fig. In the small inset of this figure the 
differences in the distribution of are shown at both sides of the transition. This distribution is peaked at = 
in the FM phase, and changes to a value clearly different from zero (around 2) when crossing the transition line to 
the MOD phase, as expected from the mean field calculation (which gives — 2 aX the transition point for all values 
of m^). 

An hysteresis might be observed for second order transitions as well when it indicates the sudden increase of the 
relaxation time around the critical point. In order to exclude this possibility we looked into a feature characteristic 
of the first order transitions only, the appearance of double-peaks in the histogram of important observables. The 
histogram of the energy, shown in the bottom part of Fig. |51 corresponds to an L = 8 simulation at C2 = 1.0, 
w? = —0.55. Such a double-peak structure is not expected when the fluctuations around the mean field solution 
of Section ImH are considered at one-loop level. It is well known that fluctuations may change the transition from 
second to first-order, especially for the so-called weak first-order transitions 0|. A first-order character is indeed 
what is expected for a liquid-solid transition where there is a finite latent heat. 

Finally, Fig. shows that the FM-MOD transition line obtained in the numerical simulation is almost horizontal, 
here in agreement with the mean field calculation, and also with the experimental phase diagram Fig. ^ 

3. MOD-PM transition 

The solid to ordinary liquid transition is also first-order in ^He. Our mean field solution predicted however a 
continuous transition. But this is again changed by the effect of fluctuations, as the numerical simulation reveals. 

The top part of Fig. 1101 shows some hysteresis plots for different lattice sizes at C2 = 1.5. Their form shows clear 
signs of a metastability at the MOD-PM transition. But notice however the appreciable finite size dependence on the 
location of the apparent transition point. This effect may be understood by looking at the shape of the distributions 

of of the L = 8 and L = 16 lattices, shown in the middle part of Fig. ^| The i = 16 distribution shows a 
peak at the value of the mean field prediction in the modulated phase (2.27 for C2 — 1.5). The L — 8 distribution 

is however extended in a larger range with a maximum at a lower value of P^. In fact, the spacing between Fourier 

modes is 27r/L, a finite number on a finite lattice, which implies that the measured values of P^ are also discretized. 
The mean field value lies between two of the allowed values for the L — 8 size, producing a competition between 
different kinds of modulation in the system. This is of course a finite size effect which disappears for large lattices. 

The P^ distribution is spread in the paramagnetic phase. This is clearly seen in Fig. llOI-middle for L = 16 while 
the spreading is weaker for L = 8. 

The bottom part of Fig. ^| shows the energy histogram of an L = 12 lattice at this transition line with the 
characteristic double-peak of a first-order transition, confirming the discontinuous character of this transition. 
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FIG. 9; Top: Hysteresis cycle in C2. Figure shows Sc, Eq. 13911 . versus C2 for — —0.55 in an L = 8 lattice. Inset: 

distribution, Eq. 1)37^ . at both sides of the transition for the same run as in the main plot. Inside the EM phase F{P^) has 

a peak in = 0, while when crossing to the MOD phase the peak changes its position to ~ 2. Bottom: Histogram of 
Sc for a run with fixed parameters C2 = 1.0, — —0.55, in an L = 8 lattice. Dotted lines mark the position of the peaks 
of the histogram in the hysteresis plot. The double-peak form of the histogram shows clearly the first order character of the 
FM-MOD transition. 

V. CONCLUSIONS 

We have proposed and studied an effective field theoretical model which is supposed to describe the ^He phase 
diagram around the point where solid, normal liquid and superfluid phases meet. This model is a generalization of 
the XY model which describes the universality class of the A-transition and it is able to explain the emergence of 
the different phases accounting for the two-excitation dispersion relation in the superfluid phase, and to relate the 
apparition of a condensate at non-zero momentum (solid phase) with a continuous deformation of the phonon-roton 
dispersion relation when increasing the pressure (here represented by the coefficients of the higher order derivatives) 
at fixed temperature from the superfluid phase. 

A mean field study, together with Monte Carlo simulations of this model, have been performed. The numerical 
simulations do not modify very much the location of the transitions in the mean field phase diagram but change in a 
qualitative way the nature of the transitions obtained in the mean field approximation. Both the form of the phase 
diagram and the order of the transitions (second-order of the XY universality class for the superfluid-normal liquid 
transition, first-order for the solid-superfluid or solid-normal liquid transitions) agree with those of *He. 

This model might be used as an starting point to study the possible apparition of other phases, such as a "supersolid" 
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FIG. 10: Top: Hysteresis cycle in m^. Figure shows Sm, Eq. 11381 . versus for C2 = 1.5 for lattice sizes of L = 8, 12, 16. 
Middle: distribution at both sides of the transition for the same run as the one in the top part for the L = 8 and L = 16 
lattices. The dotted line at — 2.27 is the value of predicted by mean field for this value of the parameter C2. Bottom: 
Histogram of Sm for a run with fixed parameters C2 = 1.5, = 0.57 in an L = 12 lattice. The histogram shows clearly the 
first order character of the MOD-PM transition. 

phase, which, after many years of debate [T^. seems to have been experimentally observed very recently [19| in the 
''He system. 
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